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Abstract 

We point out that a certain finite size effect in heavy ion physics, the 
canonical enhancement, is based on the difference of conservation con- 
strained pair statistics for Poisson and Gauss distributions, respectively. 
Consequently it should occur in a wide range of phenomena, whenever 
comparing rare and frequent events. 

It is often so, in particular in relativistic heavy ion collisions observing thou- 
sands of newly produced particles in an event, that only statistical information 
is available for tracing back characteristica of an earlier stage of interacting mat- 
ter. Aiming at the discovery of quark gluon plasma it is especially important 
to be aware of an unavoidable loss of information due to the very nature of 
statistics. Selecting out, however, relevant control parameters even statistical 
information can be used to make qualitative distinctions with respect to the 
kind of matter existed for a while in such experiments, in particular drawing 
conclusions about a possible equation of state. 

A truly interesting control parameter is the size of the reacting system, which 
is often regarded as infinite in theory (as the so called thermodynamical limit). 
In reality it is, however, finite and to a certain degree controllable in experi- 
ments by varying the target nucleus and triggering measurements by centrality 
(multiplicity, transverse energy) of a collision. But even accepting the infinite 
size assumption (theoretically in any high energy collision even infinitely many 
particle pairs can be created for a while), there is a further distinction related 
to the mean particle number in an event. 



Two famous, very characteristic limits can be considered: i) the limiting 
case of frequent events producing many particles on the mean, approximately 
described by the Gauss distribution, and ii) very rare events with few (less than 
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one) particles per event as a mean value, generally covered by the Poisson dis- 
tribution. In this letter we point out that a characteristic finite size effect in 
relativistic heavy ion collisions, named canonical enhancement, can be under- 
stood as originating in the difference between rare and frequent events, between 
sparse and cupious particle production. Mathematically this reveals itself in the 
difference between a large and a small mean particle number. The large number 
case - as it is well known - gives results equivalent to those stemming from a 
Gauss distribution (central limit theorem of statistics). 

The latter equivalence is the basis of the equivalence of the canonical and 
grand canonical approaches in the thermodynamical limit, while the Poisson 
statistics with a small mean number plots the difference between these ap- 
proaches. The ratio of the canonical and grand canonical result in small systems 
is less than one, the same result can be obtained simply relying on the Poisson 
distribution alone. The effect is in principle more general than the canoni- 
cal - grand canonical distinction; processes and phenomena featuring Poisson 
statistics occur also in physical situations far from equilibrium. Conversely, 
observation of this deviation signals truly a mesoscopic nature of the physical 
system (a fireball formed in heavy ion collisions), but does not prove in itself 
equilibrium, nor measures temperature or volume of the system. The meso- 
scopic nature reflected by the 'canonical' suppression factor can be extracted 
from mean particle numbers alone. 

The importance of this finite size effect for relativistic heavy ion collisions, in 
particular its explanation power for finding more strange particles per nucleon 
produced in ion - ion than in proton - proton collisions (the so called strangeness 
enhancement), was first realized by Redlich et. al.[[[J ^]. Originally presented as 
a constraint stemming from the proper statistical treatment of a U(l) symmetry 
in the canonical approach (whence the name canonical enhancement) ||, and re- 
ferring to volume and temperature during the derivation of the result, soon was 
it rederived on the basis of rate equations and has been shown that there can 
be an underlying master equation with the Poisson distribution as stationary 
solution. Rafelski and Letessier emphasized that this phenomenon is strongly 
related to pair statistics (to associated production of conserved charge and an- 
ticharge) and it holds also out of chemical equilibrium with general fugacity 
factors ||. Also a debate has been started about the relevance of this effect on 
the strangeness enhancement in heavy ion collisions, especially at CERN SPS 
using 40, 60 and 160 GeV/nucleon beam energies. 

Our aim in this letter is to rederive the canonical enhancement factor relying 
on the Poisson distribution alone, not assuming even thermal equilibrium in the 
background. One does not need to refer either to temperature nor to volume 
in this derivation, - even if also thermal systems in a given volume may show 
Poisson statistics, this is not their only possible origin. Then we apply this finite 
number suppression factor to an analysis of particle production in the framc- 
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work of the sudden quark coalescence picture of hadron formation, ALCOR||. 
We recover qualitatively the linear coalescence assumption for the low particle 
number case, i.e. that the number of composite objects (mesons and baryons) 
are proportional to the product of the numbers of its constituents (the quarks 
and antiquarks). 



From the probability P(n) of getting n particles of a certain kind the ex- 
pectation value (n) and higher moments like quadratic spread can easily be 
obtained by the use of the following generating function: 

Z{ 1 )=Y J P{n)e' in . (1) 

n 

The expected (mean) number becomes 

(2) 

7=0 

Under certain circumstances we are particularly interested in pairs of particles 
carrying a sum or a difference of a given physical quantity, than in the one- 
particle distribution. On the basis of independent one-particle distributions the 
pair statistics and the distribution of the sum can be calculated. Without any 
constraint on the difference, the distribution of the sum is described by the 
convolution of the corresponding distributions: 

Ms) = ^Pi(n)P 2 (m)<W„, s = ^Pi(n)P 2 (s - n). (3) 

n,m n 

This is, however, correct only if there is no information on the difference. In 
particular the production of particles and antiparticles happens always in pairs 
due to charge (and further) conservation laws, and as a consequence, even if the 
production process is statistical, their difference is bounded to be zero. In this 
case we obtain another distribution of the sum: 

This is a conditional probability with the factor Po being the probability of get- 
ting zero difference in the number statistically, which ensures the normalization 
of the result to one. Fig.l shows a geometrical interpretation of the distinc- 
tion between these two pair statistics: the unconstrained distribution can be 
obtained by adding all points in (n, to) space over a given bin of the diagonal 
representing a given sum s — n + to, while in case of the fixed zero difference 
only hits in the diagonal stripe counts (and the ratio of the stripe to the total 
area reflects the 1/Po normalization factor) . 
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Fig.l Demonstrating the constrained (main diagonal stripe) and uncostrained (the 
whole square) distribution of the sum of two Poisson deviates with mean value 10 for 
each. 



Applying this idea to restless coalescence of quarks into hadrons, we still 
constrain the difference to zero, but this time due to a confinement principle. 
The number of coalesced hadrons is the half sum in this case. Caring for pions 
only in a simplified world of light quarks and antiquarks, we arrive at having n 
from both, with zero difference and a sum of 2n forming exactly n pions. The 
generating function of the constrained distribution Pq (s|0) in this case reads as 

M7) = ^£WWe 2 ^. (5) 

n 

Assuming Poisson distribution for the quarks this generating function becomes 

Z ch) = ± E ^" W T r<9e1 " = 4-e-<«>-«> W^VWt))- (6) 
Pn n\ n\ Pq 

n 

with Iq(x) being the Bessel function of the first kind with imaginary argument. 
The normalization factor 1/Pq can be obtained from Z(0) = 1, giving 



The expectation value of the pions is 



with the label P reminding us to the Poisson distribution. For small mean 
number of quarks this leads to the product assumption of the simple linear 
coalescence model: 

(tt)p« 2(g)(5). (9) 
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For large mean numbers on the other hand the ratio of the Bessel functions, 
h/Io approaches one. This latter factor is the canonical suppression factor. 



As opposed to the above analysis of the Poisson distribution, for the Gauss 
distribution both the unconstrained and the zero difference constrained pair 
statistics leads again to a Gauss distribution (but with a larger width). As- 
suming Gauss distributed quark and antiquark numbers with respective square 
widths equal to the mean values, as it is typical for equilibrium ideal gases, the 
expectation value of the half sum, the pion number becomes the harmonic mean 
of the expectations (because by convolution of Gauss distributions the inverse 
square widths are additive): 

2(g) (g) 

(7r>G = ¥T¥' (10) 

In case of zero baryon charge (q) = (q) and these formuli reduce to a very simple 
result: 

(^g = <<?>; Wp = (9) I -r^fv v (ii) 

M 2 W» 

As a consequence one can express canonical enhancement in this case as a rela- 
tion between Poisson and Gauss expectation values. 

It is particularily interesting when we compare rare particles, for example 
K + mesons coalesced from Poisson distributed u and s quarks with copiusly 
produced n + pions glued from u and d quarks. The expected meson numbers, 

(K+) P = V^ (12) 
I (2 V us) 



for kaons and 



<* + >g = ^= (13) 
u + a 



for pions (with u, s and d denoting here the expectation values of the respective 
quark numbers) for a baryon free and strangeness free, and isotopically sym- 
metric fireball due to d = d = u = u and s = s = fd with a fixed strangeness 
ratio / = d/s reduce to 

<* + >-»^ W> (I4) 

and 

(n+) G = u. (15) 

In this simplified scenario the kaon to pion ratio can be expressed as a function 
of the pion number and so can be compared with experimental results. Of 
course the result can only be qualitative on this level, since by using the above 
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assumptions for the quark and antiquark numbers, one is bound to predict the 
same ratio for K + /ir + and for K~ /ir~ . It is not quite fulfilled in experiments: 
comparing the heavy ion results with those of proton - proton collisions, these 
ratios increase by a factor of roughly 2 and 1.6 respective lyQ. 

The result of this simple idea can be inspected in Fig. 2, where the normalized 
ratio (A' + )p/((7t + )gv / 7 is plotted versus the scaled pion number, 2(7t + )gv / 7- 
This is exactly the canonical suppression factor. Using the value / ~ 0.5 for 
estimating the Wroblinski factor / one concludes that the kaon / pion ratio falls 
to its half at an expected pion number of about 0.5 as compared to large pion 
numbers. This value corresponds to a situation before resonance decays, so it 
is not directly observable in experiments, but correcting for resonance decays is 
always possible in the framework of a tehoretical model. We used ALCOR for 
this purpose and concluded that the kaon-pion ratios normalized to the pp case 
do not change much due to resonance decay, as the unnormalized ratios do. 
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Fig. 2 Canonical suppression factor obtained from the expectation value of the sum 
of two (numerically simulated) Poisson distributions. It coincides with the analytically 
derived ratio of two Bessel functions h/Io- 



Assuming that both pions and kaons are constructed from Poisson dis- 
tributed quarks, their ratio changes between two limiting values corresponding 
to few and many quarks respectively: 

(K)p _ n iUWT) 4(2u) 

(vr)p V/ / (2«V7) h{2uY K ' 
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gives / for u -C 1 and y/J for « » 1 . Considering the light quark (pion) number 
before resonance decay the former case is realized in pp collisions (according 
to ALCOR (tt) ~ 0.1) and the latter case in PbPb collisions ((tt) ~ 10, before 
resonance decay). The double ratio, 

(K/n)p b p b _ 1 

(KM PP ~ V7 1 ; 

would be about 1.4, which is to be compared with the experimental values 1.6 
for negative and 2 for positive kaon to pion ratios. The transition is around 
(tt) w 0.7 . . . 1.4, probably occuring in a collision of relatively light ions. 

In conclusion we pointed out that the canonical suppression factor of a ratio 
of Bessel functions can be derived solely relying on properties of the Poisson dis- 
tribution, and hence is only as much related to thermal equilibrium properties 
as the latter realizes a Poisson distribution for rare particles. This distribution, 
however, may stem also from a series of dynamical events and therefore the cor- 
responding factor should show in a wide class of phenomena. We gave a rough 
estimate of reflecting this size (expected number) dependence in heavy ion col- 
lisions by the strangeness enhancement in the framework of a quark coalescence 
hadronization model, ALCOR. 
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